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We investigate the dynamics of the recently proposed model of assisted inflation. In this model an 
arbitrary number of scalar fields with exponential potentials evolve towards an inflationary scaling 
solution, even if each of the individual potentials is too steep to support inflation on its own. By 
choosing an appropriate rotation in field space we can write down explicitly the potential for the 
weighted mean field along the scaling solution and for fields orthogonal to it. This demonstrates that 
the potential has a global minimum along the scaling solution. We show that the potential close to 
this attractor in the rotated field space is analogous to a hybrid inflation model, but with the vac- 
uum energy having an exponential dependence upon a dilaton field. We present analytic solutions 
describing homogeneous and inhomogeneous perturbations about the attractor solution without re- 
sorting to slow-roll approximations. We discuss the curvature and isocurvature perturbation spectra 
produced from vacuum fluctuations during assisted inflation. 
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T. INTRODUCTION 



A single scalar field with an exponential potential is known to drive power-law inflation, where the cosmological 
scale factor grows as a cx t p with p > 1, for sufficiently flat potentials j^-Mj]. Liddle, Mazumdar and Schunck jlj] 
recently proposed a novel model of inflation driven by several scalar fields with exponential potentials. Although each 
, separate potential, 

^ Vi = V ev(-J™l+L), (1.1) 
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may be too steep to drive inflation by itself (pi < 1), the combined effect of several such fields, with total potential 
O |- energy 

n 
i=l 

leads to a power-law expansion a cx $ with Q 

b: p=5>i, (i.3) 

provided p > 1/3. Supergravity theories typically predict steep exponential potentials, but if many fields can cooperate 
to drive inflation, this may open up the possibility of obtaining inflationary solutions in such models. 

Scalar fields with exponential potentials are known to possess self-similar solutions in Friedmann- Robertson- Walker 
models either in vacuum |3|4|] or in the presence of a barotropic fluid §-§. In the presence of other matter, the scalar 
field is subject to additional friction, due to the larger expansion rate relative to the vacuum case. This means that 
a scalar field, even if it has a steep (non-inflationary) potential may still have an observable dynamical effect in a 
radiation or matter dominated era ]9L]lT|] . 

The recent paper of Liddle, Mazumdar and Schunck (jj was the first to consider the effect of additional scalar fields 
with independent exponential potentials. They considered n scalar fields in a spatially flat Friedmann-Robertson- 
Walker universe with scale factor a(t). The Lagrange density for the fields is 



- 1 

£ = E-oW - v *> (1-4) 
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with each exponential potential Vi of the form given in Eq. (LI). The cosmological expansion rate is then given by 



1 



H 2 = 



8tt 



and the individual fields obey the field equations 



(1.5) 



(1.6) 



One can then obtain a scaling solution of the form 



(1.7) 



Differentiating this expression with respect to time, and using the form of the potential given in Eq. (1.1) then implies 
that 



1 • 



= 0. 



and hence 



Pi 



The scaling solution is thus given by JTJ 



(1.8) 



(1.9) 



(1.10) 
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A numerical solution with four fields is shown in Fig. 1 as an example. In Ref. jjj the authors demonstrated the 
existence of a scaling solution for n scalar fields written in terms of a single re-scaled field = ■s/p/pifix ■ The choice 
of (pi rather than any of the other fields is arbitrary as along the scaling solution all the 0, fields are proportional to 
one another. 



FIG. 1. Evolution of four fields 
P3 = 2 and p4 = 7. 
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time, t 



and fa (from bottom to top) during assisted inflation with pi = 0.3, P2 = 1, 



In this paper we will prove that this scaling solution is the late-time attractor by choosing a redefinition of fields 
(a rotation in field space) which allows us to write down the effective potential for field variations orthogonal to 
the scaling solution and show that this potential has a global minimum along the attractor solution. In general 
the full expression for an arbitrary number of fields is rather messy so we first give, in Sect. |l[ the simplest case 



where there are just two fields, and then extend this to n fields in Sect. III. The resulting inflationary potential is 



similar to that used in models of hybrid inflation and we show in Sect. IV that assisted inflation can be interpreted 
as a form of "hybrid power-law inflation" . As in the case of power-law or hybrid inflation, one can obtain analytic 
expressions for inhomogencous linear perturbations close to the attractor trajectory without resorting to slow-roll 
type approximations. Thus we are able to give exact results for the large-scale perturbation spectra due to vacuum 
fluctuations in the fields in Sect. |y|. We discuss our results in Sect. VI. 
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II. TWO FIELD MODEL 



We will restrict our analysis initially to just two scalar fields, fa and fa, with the Lagrange density 
C 
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We define the fields 
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(2.1) 

(2.2) 
(2.3) 



to describe the evolution along and orthogonal to the scaling solution, respectively, by applying a Gram-Schmidt 
orthogonalisation procedure. 

The re-defined fields fa and a 2 are orthonormal linear combinations of the original fields fa and fa . They represent 
a rotation, and arbitrary shift of the origin, in field-space. Thus fa and <j 2 have canonical kinetic terms, and the 
Lagrange density given in Eq. (2.1) can be written as 



where 



V(a 2 ) = V 
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Pi 
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V(a 2 ) exp 
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(2.4) 



(2.5) 



It is easy to confirm that V(a 2 ) has a global minimum value Vo at a 2 — 0, which implies that a 2 = is the late time 
attractor, which coincides with the scaling solution given in Eq. Jl.lCD for two fields. 

Close to the scaling solution we can expand about the minimum, to second-order in a%, and we obtain 



V(fa,a 2 )^V 
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(2.6) 



Note that the potential for the field a 2 has the same form as in models of hybrid inflation |T!|,|2(J where the inflaton 
field rolls towards the minimum of a potential with non- vanishing potential energy density Vq . Here there is in addition 
a "dilaton" field, fa, which leads to a time-dependent potential energy density as a 2 — > 0. Assisted inflation is related 
to hybrid inflation Jl^,|2^| in the same way that extended inflation J2TJ was related to Guth's old inflation model [§2| . 
As in hybrid or ex tended inflation, we require a phase transition to bring inflation to an end. Otherwise the potential 
given by Eq. (|2.6[) leads to inflation into the indefinite future. 



III. MANY FIELD MODEL 

We will now prove that the attractor solution presented in Rcf. jl| is the global attractor for an arbitrary number of 



fields with exponential potentials of the form given in Eq. (1.1), using proof by induction. To do this, we recursively 
construct the orthonormal fields and their potential. 

Let us assume that we already have n fields fa with exponential potentials Vi of the form given in Eq. ([□]) and 
that it is possible to pick n orthonormal fields a 2 , . . . , a n and fa such that the sum of the individual potentials Vi 
can be written as 



5> 

i=l 



Pn mpi 



where we will further assume that V n = V n {a{) has a global minimum V^(0) = Vq when cr.j = for all i from 2 to n. 
It is possible to extend this form of the potential to n + 1 fields if we consider an additional field fa+i with an 



exponential potential V n+ \ of the form given in Eq. (1.1). Analogously to the two field case, we define 
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where 



Pn+1 = Pn + Pn+1 ■ 

Using these definitions we can show that the sum of the n + 1 individual potentials V, can be written as 

n+l / 



J~] Vi = V n+ i exp 



16tt 
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where V n +i = V n +i(<7i) is given by 
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(3.2) 
(3.3) 

(3.4) 
(3.5) 
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(3.6) 



Because we have assumed that V n has a global minimum value V n (0) = Vq when Oi = for all i from 2 to n, one can 
verify that V n +i also has a minimum value V^,+i (0) = Vb when a 1 , = 0, for all z from 2 to n + 1. 

However, we ha ve alr eady shown in Sect. [Il] that for two fields (fix an d </>2, we can define two fields </ >2 a nd 02, 
given in Eqs. (2.2) and (2.3) whose combined potential given in Eq. ( |2.4) is of the form required in Eq. (3.1), with 
Pi = Pi + P2- Hence we can write the potential in the form given in Eq. (3^1) for n fields, for all n > 2, with 



P=Pr, 



i=l 



(3.7) 



Equations ( |2.2|) and (3^2) then lead us to the non-recursive expression for the "weighted mean field" 



-5 # 



mpi ft , 
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(3.8) 



which describes the evolution along the scaling solution. T his i s simply a rotation in field space plus an arbitrary shift, 
chosen to preserve the form of the potential given in Eq. (3.1). The n — 1 fields ai describe the evolution orthogonal 
to the attractor trajectory. 

The pote ntia l V n h as a global minimum at Gi = 0, which demonstrates that this is the stable late-time attractor. 
From Eqs. (2.5) and (|3.6|) we get a closed expression for V n , 
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(3.9) 



Close to the attractor trajectory (to second order in <7j) we can write a Taylor expansion for the potential 

8tt 



^2Vi = V 1 



i=l 



o-j I exp 



16tt 



P TOpi 



(3.10) 



Note that this expression is dependent only upon p and not on the individual pi 
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IV. STRINGY HYBRID INFLATION 



The form of the potentials in Eqs. ( |2.6| ) and ( 3.1C| ) is reminiscent of the effective potential obtained in the Einstein 
conformal frame from Brans-Dicke type gravity theories. The appearance of the weighted mean field, (j>, as a "dilaton" 
field in the potential suggests that the matter Lagrangian might have a simpler form in a conformally related frame. 
If we work in terms of a conformally re-scaled metric 



g^u = exp 



16tt 



then the Lagrange density given in Eq. (2.4) becomes 

L = cxp 



16tt 



p mpi 



p mpi 



9p,v ■ 



n 1 



V 



(4.1) 



(4.2) 



In this conformal related frame the field <f> is non-minimally coupled to the gravitational part of the Lagrangian. 
The original field equations were derived from the full action, including the Einstein-Hilbert Lagrangian of general 
relativity, 



S 



d 4 x. 
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^R 
16tt 



(4.3) 



where R is the Ricci scalar curvature of the metric g^ v . In terms of the conformally related metric given in Eq. (4.1) 
this action becomes (up to boundary terms |f26f ) 



5' 



j4 / ~ — 4> 

a xy —ge 



n-l 



'"-PI 

16tt 



1 

R-LJ(V<S>)--Y,(VCT 1 ) 2 -V 



where we have introduced the dimensionless dilaton field 

$ = 

and the dimensionless Brans-Dicke parameter 



16tt 



p mpi 



(4.4) 



(4.5) 



(4.6) 

Thus the assisted inflation model is identical to n— 1 scalar fields cTj with a hybrid inflation type potential V(&i) in 
a string-type gravity theory with dilaton, $ oc <j>. However, we note that in order to obtain power-law inflation with 
p l the dimensionless constant u> must be much larger than that found in the low-energy limit of string theory 
where u) = — 1. 



V. PERTURBATIONS ABOUT THE ATTRACTOR 



The redefined orthonormal fields and the potential allow us to give the equations of motion for the independent 
degrees of freedom. If we consider only linear perturbations about the attractor then the energy density is independent 
of all the fields except </>, and we can solve the equation for the <f> field analytically. 

The field equation for the weighted mean field is 



+ w = x l -^^L. (5.i) 

p mpi 



Along the line <Xj = for all i in field space, we have 



V = J^JL), (5.2) 

V V P rriv\J 

and the well-known power-law solution Q with a oc i? is the late-time attractor |||J for this potential, where 

0(t) = &ln(iV (5.3) 



to 



and 4>o — m Pi\/?'/47r and to = mpi\/p/87rVb(3p — 1). 
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A. Homogeneous linear perturbations 



The field equations for the Oi fields are 



dV 
da* 



= 0. 



(5.4) 



where the potential V is given by Eqs. (3.1) and (3.9), and the attractor solution corresponds to <fi — 0. Equation ( 3.10| ) 
shows that we can neglect the back-reaction of <jj upon the energy density, and hence the cosmological expansion, to 
first-order and the field equations have the solutions 



where 



3p-l 



-1± 



' 3(p-3) 
3p-l 



(5.5) 



(5.6) 



for p > 3, confirming that ui = is indeed a local attractor. In the limit p — > oo we obtain Sj = —2. For 1 < p < 3 
the perturbations are under-damped and exe cute decaying oscillations about a% = 0. 

The form of the solutions given in Eq. (5.5) for <7;(i) close to the attractor is the same for all the orthonormal fields 
<Ji , as demonstrated in Fig. |2j. Their evolution is independent of the individual pi and determined only by the sum, 
p, as expected from the form of the potential given in Eq. ( 3.1 0| ) . 




B. Inhomogeneous Linear Perturbations 

Conventional hybrid inflation and power-law inflation are two of the very few models (l8| in which one can obtain 
exact analytic expressions for the spectra of vacuum fluctuations on all scales without resorting to a slow-roll type 
approximation. In the case of hybrid inflation, this is only possible in the limit that the inflaton field a approaches the 
minimum of its potential and we can neglect its back-reaction on the metric p3fl . As the present model is so closely 
related to both power-law and hybrid inflation models close to the attractor, it is maybe not surprising then that we 
can obtain exact expressions for the evolution of inhomogeneous linear perturbations close to the scaling solution. 

We will work in terms of the redefined fields <p and a% , and their perturbations on spatially flat hypersurfaces [p5| . 
In the limit that ct,; — > we can neglect the back-reaction of the <j{ field upon the metric and the field <j). Perturbations 
in the field </> then obey the usual equation for a single field driving inflation |l6[ , and perturbations in the field <7i 
evolve in a fixed background. Defining 

u = aS(f> , (5-7) 
Vi = addi , (5.8) 

we obtain the decoupled equations of motion for perturbations with comoving wavenumber k, 



G 



.2 



Z 



z 



u k = , (5.9) 



i& + ( fc 2 + a 2 ^ - ^ ) « tt = , (5.10) 



where |16 z = a 2 4>' /a! and a prime denotes differentiation with respect to conformal time r\ = J dt/a. For power-law 



expansion we have z oc a oc (—77) ^ and thus 



a (p-1) 

We also have ai? esc —p/((p — l)rj) which gives 

2 rf 2 V^ _ 2(3p- 1) 
° d^" ~ (p-1) 2 



a" p(2p-l) _ 2 

— T 77 . (5.11) 



rf\ (5.12) 



where we have used the fact that d 2 V/daj — IQuV/m^ along the attractor. The equations of motion therefore 
become 

<+(fc 2 - ^~ 2 1/4) )» fc = 0, (5.13) 

A 

j) 2 



< fc +(fc 2 - A2 J 1/4r )« <fc = 0, (5.14) 



where 



^3^-3)0-1/3) 
2 p — 1 

and the general solutions in terms of Hankcl functions are 

u k = Uxi-krtY'^H^i-hri) + C/ 2 (-fcr?) 1/2 ^ 2) (-^) , (5-17) 
v ik = V u (-krj) x ' 2 H£\-kTj) + V 2l {-ki 1 ) 1 ' 2 H^ ) {-ki 1 ) . (5.18) 

Taking only positive frequency modes in the initial vacuum state for \krj\ 3> 1 and normalising requires u k and 
Vik — > e~ %kv /y2k, which gives the vacuum solutions 

u k = ^(-n^e^+^H^i-kv), (5.19) 

v ik = l(-w)V* e H x+1 *H$\-kr,) % (5.20) 

In the opposite limit, i.e., — kr\ — > 0, we use the limiting form of the Hankel functions, H^{z) ~ — (i/~K)T(v)z~ v , and 
therefore on large scales, and at late times, we obtain 

u k ^ ^c^^-^'-'rw, (5.21) 

i, lfe -> ^=ie^ A (-A^)* - * T(A) . (5.22) 
V 7rfc 



The power spectrum of a Gaussian random field if> is conventionally given by = (/c 3 /27r 2 ) (It/ 1 ! 2 ). The power 
spectra on large scales for the field perturbations 5<p and 8ai are thus 
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where we have used rj — 



—{y— 1/2) /(aH) and we define 



C (a) = ^M. (5.25) 

V > 2 | r( 3 ) 

Both the weighted mean field cf> and the orthonormal fields &i are "light" fields (m 2 < 3H 2 /2) during assisted inflation 
(for p > 3) and thus we obtain a spectrum of fluctuations in all the fields on large scales. Note that in the de Sitter 
limit, p — ► oo and thus v — > 3/2 and A — > 3/2, we have T^J/ 2 — * H/2ir, and "P^ 2 — * H/2ir. 

At late times, that is fcr; — > 0, the perturbations in the weighted mean field, <50, approach a constant, while the 
perturbations in the orthonormal fields, Scfi, decay in agreement with our solutions for the homogeneous perturbations 



given by Eqs. (5.5) 



Denoting the scale dependence of the perturbation spectra by An x = dlnV x /d\nk, we obtain 

An ff * = 3-2i/ = -^-, (5.26) 
p- 1 



An { , = 3 - 2A = 3 ( 1 - VMESEUR j . ( 5 .27) 



VI. DISCUSSION 



We have shown that the recently proposed model of assisted inflation, driven by many scalar fields with steep 
exponential potentials, can be better understood by performing a rotation in field space, which allows us to re-write 
the potential as a product of a single exponential potential for a weighted mean field, (f>, and a potential V n for the 
orthogonal degrees of freedom, er^, which has a global minimum when ai = 0. This proves that the scaling solution 
found in Ref. Q| is indeed the late-time attractor. 

The particular form of the potential which we present for scalar fields minimally-coupled to the spacetime metric, 
can also be obtained via a conformal transformation of a hybrid inflation type inflationary potential jn^(| with 
a non-minimally coupled, but otherwise massless, dilaton field, oc 4>. Thus we see that assisted inflation can be 
understood as a form of power-law hybrid inflation, where the false- vacuum energy density is diluted by the evolution 
of the dilaton field. 

We have also been able to give exact solutions for inhomogeneous linear perturbations about the attractor trajectory 
in terms of our rotated fields. Perturbations in the weighted mean field <p corresponds to the perturbations in the 
density on the uniform curvature hypersurfaces, or equivalently, perturbations in the curvature of constant density 
hypersurfaces: 

C=^. (6.1) 

4> 

These perturbations are along the attractor trajectory, and hence describe adiabat ic cu rvature perturbations. The 
spectral index of the curvature perturbations on large scales is thus given from Eq. ( 5.26| ) as 



and is always negatively tilted with respect to the Harrison-Zel'dovich spectrum where n s = 1. Note that in the 
slow-roll limit (p — > oo) we recover the result of Ref. jj]. 

First-order perturbations in the fields orthogonal to the weighted mean field are isocurvature perturbations during 
inflation. Vacuum fluctuations lead to a positively tilted spectrum. The presence of non-adiabatic perturbations 
can lead to more com plic ated evolution of the large-scale curvature perturbation than may be assumed in single- 
field inflation models |27|-|29(| . However, we have shown that these perturbations decay relative to the adiabatic 
perturbations and hence we recover the single field limit at late times. In particular we find that the curvature 
perturbation £ becomes constant on super-horizon scales during inflation. Note, however, that assisted inflation must 
be ended by a phase transition whose properties are not specified in the model. If this phase transition is sensitive to 
the isocurvature (non-adiabatic) fluctuations orthogonal to the attractor trajectory, then the curvature perturbation, 
C, during the subsequent radiation dominated era may not be simply related to the curvature perturbation during 
inflation. 
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